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Most exoplanetary systems in binary stars are of S-type, and consist of one or more planets 
orbiting a primary star with a wide binary stellar companion. Gravitational forcing of a 
single planet by a sufficiently inclined binary orbit can induce large amplitude oscillations of 
the planet’s eccentricity and inclination through the Kozai-Lidov (KL) instability |[liS]* KL 
cycling was invoked to e^lain: the large eccentricities of planetary orbits [0]; the family of 
close-in hot Jupiters and the retrograde planetary orbits in eccentric binary systems 
10,0]* However, several kinds of perturbations can quench the KL instabilitv, by inducing 
fast periapse precessions which stabilize circular orbits of all inclinations yt): these could 
be a Jupiter-mass planet, a massive remnant disc or general relativistic precession. Indeed, 
mutual gravitational perturbations in multiplanet S-type systems can be strong enough to 
lend a certain dynamical rigidity to their orbital planes [g]. Here we present a new and faster 
process that is driven by this very agent inhibiting KL cycling. Planetary perturbations 
enable secular oscillations o^lanetary eccentricities and inclinations, also called Laplace- 
Lagrange (LL) eigenmodes 1191 ]. Interactions with a remnant disc of planetesimals can make 
planets migrate, causing a drift of LL mode periods which can bring one or more LL modes 
into resonance with binary orbital motion. The results can be dramatic, ranging from excita¬ 
tion of large eccentricities and mutual inclinations to total disruption. Not requiring special 
physical or initial conditions, binary resonant driving is generic and could have profoundly 
altered the architecture of many S-type multiplanet systems. It can also weaken the multi¬ 
planet occurrence rate in wide binaries, and affect planet formation in close binaries. 


The fiducial system has two planets on initially coplanar orbits around a solar mass primary star: 
an interior lOMjup planet on a circular orbit with an initial semi-major axis of al^ = 5 AU, and 
an exterior lOM^ planet with initial semi-major axis between 8 and 11 AU and eccentric- 
The binary is also a solar mass star with semi-major axis ai, > 100 AU and 


ity Ct = 0-05. 


corresponding period Tf, and angular frequency rib. Planetary migration driven by scatterin 


of 

planetesimals has a long and productive history in relation to solar system archeology [Il0l - ll2il . It 
is a complex process, as discussed in the Supplementary Notes. Here we use it in a simple man¬ 
ner: the outer planet is allowed to migrate outward due to interactions with a planetesimal disc, 
with its semi-major axis having a prescribed form, with characteristic timescale r. For the solar 
system, there are plausible arguments that the migration time r > 10® yr iHHi . with a lower bound 
r > 10^ yr arguably needed to recover the properties of the Neptune Trojans 11131] : we assume 
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2.5 X 10^ yr < r < 5 X 10® yr. The physics of the problem appears clearest in a secular setting 
wherein the fast planetary (but not the binary) orbital motions are averaged over, turning a point 
mass planet into a shape and orientation changing Gaussian wire We present a numerical 
simulation with a state-of-the-art, iV-wire algorithm [Il4ll . and develop a mathematical model to 
understand the results. 


In the A^-wire experiment the binary orbit was coplanar and circular, with = 1000 AU, 
and a period Tb = 22.36 Kyr. The initial periods of the two LL eigenmodes are 3.53 Myr (a 
slow mode determined mainly by the massive inner planet) and 12.25 Kyr (a fast mode reflecting 
mainly the precession of the outer planet). Outward migration of the outer planet slows down the 
faster LL mode until its period approaches the binary period Tb = 22.36 Kyr. Fig.la shows that 
the eccentricity of the outer planet eout — 0.05 until Oout — 11-89 AU. Then it is captured into 
a resonance and Cout begins increasing, rising to 0.51 when Oout — 15 AU at time t ~ 4055 Tb. 
Capture is also apparent in the behavior of the resonant argument, 0res(f) = ^out(f) — ribt, where 
Wont is the apsidal longitude. From Fig. lb we see that, after a period of circulation, 0res enters into 
libration at resonance passage, with libration maintained for the full duration of the simulation. 
More details of the capture process are given in Extended Data Figs.l(a,b). 


This capture phenomenon is closely related to the lunar evection resonance, that may have 
played a significant role in shaping the early history of the lunar orbit [ 15. H- However, the LL- 
mode evection resonance (LLER) is a new process, so we also present an analytical model, valid 
for arbitrary binary eccentricity, in the Supplementary Notes. Eor a circular binary orbit, EEER 
dynamics is governed by the normal form Hamiltonian of eqn(l25l): 




+ rf" 
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+ rf 
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where rj and ^ are a canonically conjugate pair of LE mode variables, and the parameters 6, a' 
and /S' are functions of the slowly migrating Oout- The theoretical prediction for the location of 
the exact resonance is shown as the dashed red curve in Eig.la: exact resonance is first met by a 
zero eccentricity planet around Oout — 11-875 AU; the planet circulating at eout — 0-05 encounters 
resonance a bit later (when Oout — 11-89 AU), then gets engulfed by a growing and migrating 
nonlinearly bounded resonance region. Prediction follows simulation in the mean until Cout — 
0.2 in our 4th order model; a higher-order expansion will improve the fit between model and 
simulation. The evolving topology of flows in the phase space, along with key structural 
features in and around resonance are discussed in the Supplementary Notes. 


Whereas encounter with the EEER is certain with migration, capture in it is probabilistic, and 
depends on the strength of the resonance, the migration rate and the initial planetary eccentricity 
at which EEER is encountered. Eor r ~ 10"^ T;, and initial Cout = 0.05, the probability of capture 
in EEER exceeds one-half. At the assumed binary separation, capture becomes certain as the 
migration rate is slowed down by two orders of magnitude; more details are discussed in the 
Supplementary Notes. Capture is likely to improve in tighter S-type systems: for Ob = 200 AU 
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and Oin = 5 AU, evection is crossed at Oout — 7.43 AU, with faster apse-preeession and a stronger 
resonanee, in the eourse of migration. The eapture probability remains high for initial Cout = 0.05, 
even at faster migration rates. Having deseribed the broad seeular skeleton of LLER, we note that 
the full problem is rieher due to the interplay of planetary mean-motion resonanees (PMMR). To 
study this, it is neeessary to perform simulations that do not average over planetary mean-motions. 
Below we present two sueh A^-body simulations with the open-souree paekage MERCURY lll2 . 


m 


The first MERCURY experiment is the unaveraged version of the A^-wire simulation of 
Eig.l, and its results are displayed in Eig.2. Signs of PMMR are apparent in Pig.2a, in the jumps 
experieneed by semi-major axis of the outer planet as it migrates. The system is eaptured in 
PEER with eonsequent growth of the eeeentrieity (Pig.2b), and libration of the resonant argument 
(Extended Data Pig.2a). What is remarkable though, and distinet from the A^-wire experiment, 
is the non-monotonie behavior of the mean eeeentrieity of the EEER-loeked planet, leading to 
eseape from PEER altogether, eventually settling at Cout — 0.12. Eseape from eapture is due to 
planetary mean motion resonanees (PMMR), whieh enhanee exehange of angular momentum. Of 
the four PMMR loeated near Oout — 12.56 AU the strongest is the 4:1, with argument 04:i = 
4Aout — Ain — SccJout; a short time segment of <p 4 -,i is shown in Extended Data Pig.2b. 

In the seeond MERCURY experiment, the binary orbit had eoeentrieity 0.4 and inelination 
40°, whieh are modest values for wide-binaries. Pig.3a shows a 3:1 PMMR exeiting Cout to 0.1 at 
t ~ 21 Myr, with eapture in PEER at f ~ 72 Myr when Oout — 11.92 AU. As earlier, passage 
through the 4:1 resonanee forees the outer planet out of PEER. This is followed by another exeita- 
tion around 120 Myr, assoeiated with passage through a 9:2 resonanee; and then through a eluster 
of resonanees around 14.63 AU. Meanwhile Oout grows with jumps at the PMMR (Extended Data 
Pig.3a), and 0res transits in and out of libration during PEER (Extended Data Pig.3b). Both Oout 
and Cout diffuse until the planet is ejected from the system altogether. Ejeetion is not a neeessary 
outeome, but is often assoeiated with PMMR when both Oout and Cout are large. In Pig.3b we follow 
the exeitation of the mutual inelination to 12°, due to eoupling within EEER-loek, of eoeentrieity 
and inelination by a vertieal resonanee, whieh is followed by another exeitation to 14.7°. 


PEER is a powerful and generio meohanism that oan profoundly affeot the arohiteoture of 
multiplanet S-type binary systems. It oan also oome in different flavors. 1. Inward migration 
of the inner planet oan ooour through a runaway prooess H, whose slower migration rate has 
higher eapture probability, partioularly in tighter S-type systems with shorter EE periods. Inward 
migration may explain the largish eooentrioities and inolinations in systems with super-Jupiter 
sized planets on sub-AU orbits. 2. EEER-indueed disruption in moderately wide binaries (ob < 
1000 AU) may be responsible for the reoently reported dearth of multiplanet systems in binaries 
at sueh separations lll9l] . The extent to whieh PEER disrupts/suppresses planet formation when 
Oft < 20 AU ll^ needs to assessed within planet formation studies [21]. 3. A multi-mass planetary 
system will have a broader speetrum of EE frequeneies than a two-planet system. The rieher PEER 
and stronger PMMR open more pathways for disruption, and eould relieve an initial multiplanet 
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system of all but one of its planets. 
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Figure 1: Capture into the Laplace-Lagrange Evection Resonance. The fidueial iV-wire ex¬ 
periment was performed with exponential migration, aout(^) = flout exp[f/r], with 0 *^^ = 10 AU 
and r = 10^ T;,. (a) Growth of Cout when it is eaptured in the migrating LLER. The dashed red line 
is the predietion from the analytieal 4th-order theory presented in Supplementary Notes, (b) 
transitions from eireulation to libration around 90° when eaptured in LLER. 















Figure 2 
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Figure 2: LLER with PMMR for coplanar circular binary orbit. The first MERCURY experi¬ 
ment was eondueted with damped migration, aout(f) = Ofin — Aaexp[—t/r], with Ofin = 15 AU, 
Aa = 3 AU and r = 10® yr = 4500 T^. (a) Migration of Oout with signs of PMMR. (b) LLER is 
eneountered around Oout — 11-88 AU at f ~ 15 Myr, with initial growth of Cout during eapture, 
then deeay beeause of interruption by PMMR, and ultimately eseape from LLER. 
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Figure 3 
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Figure 3: LLER with PMMR for inclined and eccentric binary orbit. The second MERCURY 
experiment was conducted with exponential migration aout(i) = ^out exp[t/r], with 0 *^^ = 10 AU 
and r = 4.4 x 10® yr = 2 x lO^T;,. (a) The 3:1 PMMR around 21 Myr excites Cout to about 0.1. 
LLER is encountered around Oout — 11-92 AU, at t ~ 72 Myr. (b) Mutual inclination first excited 
to 12°, then to 14.7°. 
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Supplementary Notes 


A Simple Analytical Model of the Laplace-Lagrange Evection Resonance 

We present an analytieal model of the Laplaee-Lagrange Eveetion Resonanee (LEER) in a two- 
planet system orbiting the primary star, with the eompanion star orbiting in the same plane as 
the planets. The main general result is eqn. dlTT) for the seeular Hamiltonian for the two Eaplace- 
Eagrange (EE) modes (given in the mode aetion-angle variables to 4th order in the planetary ee- 
eentrieities) when they are foreed by a binary orbit of arbitrary eceentrieity. In order to understand 
the fidueial iV-wire simulation reported in the main text we speeialize to a eireular binary orbit. 
Then the seeular Hamiltonian reduees to the normal form, of eqn. (l25]) . Simple eomputations 
with i/nf provide (a) a graphie narrative of the unfolding of EEER in phase space (Extended Data 
Eig.4); (b) characteristics of the EEER islands including measures of the adiabaticity (Extended 
Data Eig.5); (c) the dependence of capture probability on the initial planetary eccentricity and 
non-adiabaticity. 


The Hamiltonian governing the secular dynamics of planets of mass mi and m 2 can be 
written as: 


H,. 


—Gmim2 



GMettii 


1 

|n(f) - ri 


rb{t)-ri\ 

nity / 


— G ]\dQTn2 


1 

\rb{t) - r2 


rb{t)-r2\ 

nitr / ■ 


( 1 ) 


Here means that the expression inside is to be averaged over the Kepler orbits of both 

planets, and “< >” means that the averaging is to be performed over the Kepler orbit of either 
planet 1 or 2, as the case may be. Eet (oi, 02, at) and (ei, 62, Cb) be the semi-major axes and the 
eccentricities of planet 1, planet 2 and the binary orbit, respectively. Eet gi and g 2 be the periapse 
angles of the orbits of planets 1 and 2 , and 9 b{t) be the polar angle to the location of the binary 
star. In the absence of planetary migration, secular dynamics conserves all the semi-major axes. Cb 
is constant, because the binary is assumed to be in a given Kepler orbit. The secular Hamiltonian 
governs the dynamics of the quantities {gi, g 2 ; ei, 62 ). We assume that Oi < 02 ab, and expand 
the binary potential to quadrupolar order: for z = 1 , 2 , 


\rb{t) - Vil 


rbit)-ri\ 

nitr / 


,2 r 


4r^ 


. 3 2 15 

1 + + y (2^* 



+ ... . 


( 2 ) 


However, a multipolar expansion of the interaction between planets 1 and 2 may not be appropriate, 
because their semi-major axes may be of comparable magnitudes. Therefore we assume that both 


9 











ei -C 1 and 62 1 and expand to fourth order in the eeeentrieities: 

CLi 102 1 / \ 02 04 

= ^1^2061 + Cii 6162 cos (^1 - ^2) + Co2e2 + 0406^ 

+ C34 6462 cos {gi - g2) + C22 6462 + C22 6462 cos { 2 gi - 2 g 2 ) 

+ cjg 646^ cos {gi - g2) + + ... | , ( 3 ) 

where the c™; are funetions of ( 04 / 02 ), and ean be written in terms of Laplace coefficients [@]. 
When eqns. dH) and O]) are substituted in eqn. ([T]), we obtain the secular Hamiltonian as a function 
of the four dynamical quantities {gi,g 2 ; 61 , 62 ). However, these quantities are not canonically 
conjugate variables. Therefore we define new canonical coordinates (gi, 52 ), and their canonically 
conjugate momenta (pi,P 2 ) by, 

Qi = - \l2mi (GM^Oi)^^^ [1 - (1 - e 2 )V 2 ] sin g ^, 

Pi = + ^2mi (GM^Oi)^/^ [1 - (1 - e-)i/2] cos^i. (4) 

Expressing Q and dU) in terms of the variables {qi,Pi), the secular Hamiltonian of eqn. dD can be 
written as the sum of three terms: 





LL 


ifl 


bin 


Hr. 


(5) 


Here is the Laplace-Lagrange Hamiltonian that consists of all the time-independent quadratic 
terms, if bin is the time-dependent driving due to the binary motion, and iinon is the nonlinear part 
that has all the time-independent fourth order terms. 

iiLL = ai{ql+pl) + HQiH + P1P2) + 02(52+^2) > ( 6 ) 


where the coefficients ai, (3 and 02 are defined by. 
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fibL determines the two LL modes of oscillations of the eccentricities and periapses of the two 
planets. It has contributions from planetary interactions and the orbit-averaged binary quadrupole. 
To quadratic order the purely time-dependent binary forcing is represented by: 


Tibin = a6i(f) [ql +pV\ + 761 (b {2giPi sin [26'b(t)] + (g? - p?) cos [2di,{t)]] 

+ ab 2 {t) [ql +P 2 ] + 762 (f) { 2 g 2 P 2 sin [26>b(f)] + {ql - pt) cos [26'b(f)]} , (8) 
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where the eonstants abi(t), 7 bi(t), ab 2 {t) and 762 (^) are defined by, 
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( 9 ) 


When the binary orbit is eireular, and = 0; then the eoeffieients abi and ^62 both vanish. 

The nonlinear, time-independent, fourth-order nonlinear terms are gathered together in: 




Vi {qI +pIY + (9i92 + P 1 P 2 ) {ql +pD + P {ql + pI) {ql + pI) 

k -2 {qiq 2 + P1P2) [ql + pI) + A {qiq 2 + pmf + P 2 [ql + plf , (10) 


where the eonstants 71 , ki, p, ^ 2 , A and 72 are defined by. 
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Hnon determines the response of the LL modes to the resonant foreing by the binary. The see- 
ular Hamiltonian, defined by eqns. ([5])— (fTTl) . governs the dynamics of the LL modes, in the 4- 
dimensional phase space, (gi, ^2 ;Pi,P 2 )- Below we write each of its three terms, LTll, -^bin and 
LTnonj in terms of the action-angle variables of the LL modes. 


la. Laplace-Lagrange modes: Define new canonical variables, (Qi, Q2; -Pi, -P2) by: 

Ql = qicosx + q2smx, Q2 = - gi sinx + ^2 cosx ; 

Pi = picosx + P2sinx, P2 = -pi sin X + P2 cos X, ( 12 ) 

where x is such that tan (2x) = /3/{ai — 02 ). In the new variables, Pll of eqn ® is: 

-Pll = [tti cos^X + /^ sinXcosX + ^2 sin^x] {Qi + Pi) 

+ [ai sin^ X - /? sinXcosX + 02 cos^ x] {Ql + P2) ■ ( 13 ) 


11 


























We now define aetion-angle variables for the LL modes, (Ji, J 2 ; ^’ 2)5 by: 

Qi = sin - 01 , Pi = 008 01 , Q2 = sin 02, P2 = cos02- ( 14 ) 

Then 

-f^LL = CJi Ji + UJ2J2, (15) 

is in eanonieal form where 

CJI = (01 + 0:2) + (oi - 02) cos (2x) + / 9 sin( 2 x), 

OJ2 = (01 + 02) - (oi - 02) cos (2x) - /Ssin(2x), ( 16 ) 

are the mode frequeneies. 


lb. Resonant driving: We use eqns. (O and (fT4l) to work out ifbin in terms of aetion-angle 
variables for the LL modes. Dropping the fourth order terms in eqn. ([8]), we have 

Hun = +l(f) Jl + +2(1) J2 + +3(f) \/JlJ 2 COS (01 - 02) + Vl{t) Ji COS [20i + 2001)] 

+ 1^2(1) J2 COS [202 + 29 b(t))] + 1/3(1) J1J2 cos [01 + 02 + 29 b(t)] , ( 17 ) 

where the new eoeffieients, /ii(l), +2(1), +3(1), ^^i(l), 1^2(1) and 1/3(1), are defined by 

+i(f) = + 2 ofei(l) cos^x + 2062(1) sin^ X , 

+2(f) = +2obi(l) sin^x + 2 o 62(1) cos^x, 

/i3(l) = 2 [062(1) - 06i(l)] sin ( 2 x), 

1/1(1) = - 2761 ( 1 ) cos^x - 2762(1) sin^X, 

1/2(1) = - 2761 ( 1 ) sin^x - 2762(1) cos^X, 

1/3(1) = -2 [762(1) - 7w(i)] sin(2x) • ( 18 ) 


Ic. Secular nonlinearities: Lastly, we write H^on in terms of LL-modal variables by using 
eqns. (fT^ and (fT4l) to susbtitute for (gi,g2 ;Pi,P2) in terms of (0i,02 ; Ji, J 2 ) in eqn. (fTO . Of 
the many terms, those proportional to cos (0i — 02) and cos (20i — 202) can be dropped when 
cji and a;2 are well-separated (i.e. non-degenerate, as in the example explored in the body of the 
article), because the angle-dependent terms are oscillatory and do not contribute significantly to 
the dynamics. Therefore, nonlinear part of the Hamiltonian for nondegenerate modes can be taken 
as, 

= ilJl + i2Jl + izJlJ2, (19) 
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where the eoeffieients, ^2 and ^3 are given by 


= dr^icos'^x + 2 ki cos^ X sin(2x) + psin^(2x) 

+ 2^2 sin^ Xsin(2x) + Asin^(2x) + 4 ? 72 sin^X) 

^2 = 4 r 7 isin"^x ~ 2 ki sin^ x sin( 2 x) + psin^( 2 x) 

— 2^2 cos^ Xsin(2x) + Asin^(2x) + 4r72Cos'^X) 

.^3 = 4 ? 7 isin^( 2 x) — 2Kisin(4x) + 4pcos^(2x) 

+ 2K2sin(4x) + 2Acos(4x) + 47^2 sin^(2x) • (20) 


2. Secular Hamiltonian for nondegenerate LL modes with binary driving: Gathering together 
with the expressions in eqns. (fTSl) . (ITT]) and (fT9l) . we have the seeular Hamiltonian in the desired 
mode variables: 


Hsec = [^1 + Jl + [uJ2 + fj'2(t)] J2 + + '^2 <^2 + 

+ v J1J2 cos {'i/ji - - 02 ) + Ji COS [ 2 ^jJl + 29 h{t)] + 

+ l^2{t) J2 COS [ 2'^2 + 20 ft(t))] + y/ J1J2 cos ["^i + "^2 + ■ ( 21 ) 

There are resonanees between the binary and the LL modes, when is eommensurate with any of 
the frequeneies Ui, U 2 or (wi + ci;2)/2. The set of resonances is particularly rich for an eccentric 
binary orbit. When the binary orbit is circular, the coefficients Hiit) all vanish, and the i>i become 
time-independent. Then eqn. (l2TI) simplifies to: 

Hcivc = cci Ji + U2J2 + iiJi + ^2 >^2 + <^1 cos ( 2 '^i + 2 nfet) 

+ 1^20 J2 cos (27/>2 + 2 nbt) + 1/30 V J1J2 cos + 7/^2 + 2 n 6 f), ( 22 ) 


where the coefficients, i^io, 1^20 and z/30, are given by 
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tti sm X + 0,2 cos X 


sin(2x). 


(23) 


It turns out that oji and U 2 are negative, giving rise to three types of LLER, for rib — |a;i|, or 
rib — \^ 2 l or rib \ (<^i + a;2)/2| . The possibilities are extremely rich, so we focus on the case 
relevant to the fiducial iV-wire experiment described in the main text. 
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3. Normal form Hamiltonian for the N-wire experiment: In the A^-wire experiment, the inner 
planet has mi = 10 M^^p, Oi = 5 AU with initial ei = 0, and gi = 0; the outer planet has 
m2 = 10 M®, 02 < 11.89 AU (here it is initially at 11 AU) with 62 = 0.05 and (72 = 0. The binary 
eompanion is also a solar mass star, on a eireular orbit with semi-major axis at, = 1000 AU and 
period Tb = 22, 360.69 AU. The two LL mode frequeneies are cui = —1.33 x 10“® rad/yr (a period 
of 4.73 Myrs) and UJ 2 = —3.95 x 10“^ rad/yr (and a period of 15, 918.29 yr). Here we study LLER 
when Ub ~ |a;2|. Sinee cji and 002 are well-separated in magnitude, it is elear that rib cannot be 
elose to either |a;i| or | (cji + a;2)/2|. Then the driving terms proportional to cos (2'0i + 2nbt) and 
cos (^1 + V’2 + 2nbt) are oscillatory, and ean be dropped. Hence ETcirc is effectively independent 
of the angle which implies that Ji = Jio — constant. Therefore the resonant Hamiltonian for 
describing LLER of the second mode takes the simple form, 

Hm2 = {^2 -\-nb + ^2, Jio) J 2 -\- ^2 J 2 + ^20 J 2 COS ( 2 (^ 2 ) • ( 24 ) 


In the absence of planetary migration, this is a time-independent 1 degree-of-freedom Hamilto¬ 
nian in the canonically conjugate variables, <^2 = //^2 + ^bt and J2, and the dynamics is obviously 
integrable. This Hamiltonian is typical of 2nd-order resonance models, and can be further reduced 
to a normal form in the new canonical variables, ^ = \/ 2 J 2 cos(<yC2) and g = \/2J2 sin(<y52): 


where 




6 — (cu 2 -f rift-f ,^3 Jio) , 


a ' = -^2, 


/ 3 ' — — 1 ^ 20 - 


(25) 

(26) 


The normal form Hamiltonian has a long history in solar-system dynamics (see ll22l] and references 
therein). Of relevance to our problem is lunar evection, the resonance between the precession of 
the peripase of the Moon’s orbit around an oblate Earth, and the mean motion of a massive outer 
perturber, the Sun lllSn . 


4. Planetary migration: Numerical simulations of planetary migration with gaseous discs have 
reported a wide variety of behaviour — planets opening gaps, clearing out inner discs, stalling in 
their migration, reversing migration; multiple planets undergoing divergent migrati on, undergoing 
convergent migration, or getting captured into mean motion resonances ll23l - l25n — but we do 
not explore this here. We consider planetary migration driven by scattering of planetesimals. 
This process is believed to have taken place in the solar system, and to have left its signature 
in the dynamical properties of minor bodies, as well as spin and orbital features of the planets 
themselves lliol ll iL 126142811 . The planetary system is considered fresh out of the evaporation of 
the gaseous disc, with a remnant disc of surviving planetesimals which, in the course of their 
dynamical stirring, then scattering, by the planets is expected to drive migration. The timescale 
of migration is set by the inner boundary, mass and size distribution in the remnant disc 1111,1^. 
Eor the solar system, there are plausible arguments for it being on the order of a few times 10® 
years [Q , with a lower bound of a few 10^ yr set by exercises which seek to recover properties 
of Neptune Trojans with planetary migration Q. We cannot commit to any particular timescale 
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without careful (and largely numerieal) treatment of mean motion resonanees and their stirring of 
a preexisting dise into planet erossing orbits. We have thus assumed a range of timeseales 10^ yr 
to 4 X 10® yr, whieh is about 500 to 20, 000 binary orbital periods. The direetion of migration on 
the other hand is largely determined by the mass and loeation of the perturbing planets. 

5. Dynamics of LLER: As 02 inereases from an initial value of 11 AU due to planetary migration, 
the eontrol parameters (5, a',/3') aequire slow time dependenee, making of eqn.(l251) a 1.5 
degree-of-freedom system. In the A^-wire fidueial system, 5 is an inereasing funetion of time, 
starting with a negative value —1.14 x lO”"^ and then transitioning to positive values around 02 = 
11.884 AU; f5' ~ 10“® is always positive; a' = 2.24 initially, and deereases while remaining 
positive over the relevant range of 02. The variation of (S, a\ (3') results in non trivial ehanges 
in the topology of the instantaneous global phase portraits of i/nf in the plane. As ean be 
seen in the four panels of Extended Data Fig.4, the origin (0,0) — eorresponding to a eireular 
orbit — is always an equilibrium point. Sinee a' > 0, the origin is initially stable beeause 5 < 
—/S'. As 5 inereases in the eourse of migration, it goes unstable for 5 > —{3', whieh happens 
at 02 = 11.875 AU. This first bifureation gives rise to two stable equilibria at (ibr^cO), where 
Vc = [(<5 + l3')/<y']^^‘^. These are the eentres of LLER with librating orbits around them; see 
Figs.SI (a, b). Post-eneounter, as 6 eontinues to inerease, the eentres drift apart and the islands 
grow, eapturing into LLER any trajeetory that eomes their way. Stability is restored to the origin 
for 6 > +f3' (at 02 = 11.89 AU). This seeond bifureation gives rise to two unstable equilibria 
at (0,±^un) where ^un = [(<5 —/9')/a']^'^^; see Extended Data Figs.4(e, d). As 6 eontinues to 
grow, the basin of eireulating orbits around the origin also grows, squeezing the LLER islands and 
eapturing some of their librating orbits. Extended Data Fig.5a shows the evolution of r]c, and the 
extrema of the separatrix. In Extended Data Figs.5(b,e) we map the evolution of rjc to that of the 
planetary eeeentrieities. The dashed red curve in Fig. la of the main text is obtained from € 2 ( 3 ) of 
Extended Data Fig.5b. 


When (5, a', ( 3 ') vary slowly eompared to the libration period around LLER, we are in the 
adiabatic regime. At any time, there is a maximum eeeentrieity, emax(^)> that is reaehed by the 
separatrix; let Cc be the maximum value of all the eraax{t)- Capture is eertain if LLER is eneoun- 
tered when 62 < Cc. In our problem, Cc — 0.054, eorresponding to the the onset of the seeond 
bifureation shown in Extended Data Fig.4e. If 62 > Cc at resonanee eneounter, eapture is not 
eertain. The probability of eapture ean be eomputed analytieally ll22ll . and is given by the ratio 
of (a) the rate of inerease of the area of the libration zone, to (b) the rate of inerease of the sum 
of the areas of the libration and eireulation zones. Note that the eireulation zone has zero area 
for —( 3 '/a' < 5 /a' < ( 3 '/a', henee eapture is eertain with 5 /a' inereasing past —ft'jo/! . When 
the variation is non-adiabatic, outeomes are not easily predietable from the instantaneous phase 
portraits. Then eapture and eseape must be quantified through eomputations with H^f for different 
initial eonditions of the planet. 


6. Estimates of Adiabaticity and Capture Probability: We eomputed T = a;iib/27rrmig, where 
= d\nrjc/dt is the migration rate of the island eentres, and ccub = 2 ^y 13' [S + 13'] is the 




libration frequency around the island centre. F is a measure of adiabaticity, and is plotted versus 
Bmax in Extended Data Fig.S.d. The dynamics is increasingly adiabatic for larger 02, with the 
island migrating to higher eccentricities. For r ~ 10"^ Tf, and Cmax = 0.054, we have F ~ 0.2; 
the migration rate is larger than the libration time, implying non-adiabatic passage (in the range 
of Bmax for which capture is guaranteed in the adiabatic limit). Were r larger by a factor 100, we 
would be in the adiabatic regime, and capture in FFER would be certain for 62 < 0.054 (excepting 
near-zero eccentricity where adiabaticity is practically impossible). Theoretical estimates of the 
capture probabilities are not well-determined in this non-adiabatic regime. Hence we integrated 
trajectories with the evolving for a range of initial eccentricities and uniformly distributed 
periapses, and discovered that: (a) Capture is ruled out for 62 < 0.028; this outcome appears 
typical of non-adiabatic passage through FFER, and was already noted in studies of the early 
history of the lunar orbit [Q; (b) Matters improve for larger eccentricities: more than half the 
planets with 62 = 0.05 get captured, and the capture probability rate gets closer to the adiabatic 
estimate with larger 62 at encounter. 


References 

22. Borderies, N. & Goldreich, R A simple derivation of capture probabilities for the J + 1 : J 
and J + 2 : J orbit-orbit resonance problems. Celestial Mechanics 32, 127-136 (Feb. 1984). 

23. Alexander, R. D. & Armitage, P. J. Giant Planet Migration, Disk Evolution, and the Origin 
of Transitional Disks. Ai'trop/zy^'. J. 704,989-1001 (Oct. 2009). 

24. Armitage, P. J. Dynamics of Protoplanetary Disks. Anna. Rev. Astron. Astrophys. 49, 195- 
236 (Sept. 2011). 

25. Hasegawa, Y. & Ida, S. Do Giant Planets Survive Type II Migration? Astrophys. J. 774, 146 
(Sept. 2013). 

26. Fernandez, J. A. & Ip, W.-H. Some dynamical aspects of the accretion of Uranus and Neptune 
- The exchange of orbital angular momentum with planetesimals. Icarus 58, 109-120 (Apr. 
1984). 

27. Hahn, J. M. & Malhotra, R. Orbital Evolution of Planets Embedded in a Planetesimal Disk. 
Astron. J. 117, 3041-3053 (June 1999). 

28. Morbidelli, A., Fevison, H. F, Tsiganis, K. & Gomes, R. Chaotic capture of Jupiter’s Trojan 
asteroids in the early Solar System. Nature. 435, 462-465 (May 2005). 


16 


Extended Data Figure 1 
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Extended Data Figure 1: Capture into the Laplace-Lagrange Evection Resonance. The fidu¬ 
cial A^-wire experiment was performed with exponential migration, aout{t) = exp[l/r], with 
= 10 AU and r = 10"^ T^. (a) Phase space trajectory of the outer planet during capture. 
r] = [1 — \/l — CouJ sin(0i.es) and ^ = [1 — y/T^^^Iut] cos(0res] are canonical coordinate and 
momentum of the captured EE mode. Trajectory clearly reveals transition from circulation with 
initial eccentricity of 0.05 (inner ring) to libration when captured in EEER (funnel from inner ring 
moving toward negative t]). (b) Capture is also seen, albeit in a subdued manner, in the modest 
growth of ein during EEER. 












Extended Data Figure 2 
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Extended Data Figure 2: LLER with PMMR for coplanar circular binary orbit. The first 
MERCURY experiment was eondueted with damped migration, aout(f) = Ofin — Aaexp[—f/r], 
with Ofin = 15 AU, Aa = 3 AU and r = 10® yr = 4500 Ti,. (a) Transitions of 0res from eireulation 
to libration during eapture and baek to eireulation after eseape. (b) A 100 Kyr time-segment of 
04:1 when still eaptured in eveetion; signature of the PMMR is evident in the repeated transitions 
between libration and circulation. 
















































































Extended Data Figure 3 
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Extended Data Figure 3: LLER with PMMR for inclined and eccentric binary orbit.The second 
MERCURY experiment was conducted with exponential migration aout(i) = ^out 
®out = 10 and r = 4.4 x 10® yr = 2 x lO^Tb. (a) Migration of Oout with signs of PMMR 
around 21, 90,120 and 125 Myr. (b) Transitions of cp^es from circulation to libration during capture 
in EEER, and then back to circulation due to passage through the 4:1 PMMR which forces the outer 
planet out of EEER. 




























































Extended Data Figure 4 
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Extended Data Figure 4: Phase Space with Migrating Planet. Isoeontours of at different 
times, showing bifureations of equilibria and emergenee in islands where eapture is probable. 
Note: both ^ and r] have been resealed by a faetor of \/m 2 \/GMAa 2 to turn them into eoeentrieity 
like variables, (a) At 02 = 11 AU the origin is stable with eireulating orbits around it. (b) 
At 02 = 11.88 AU the origin has gone unstable, and two EEER islands have appeared, (c) At 
02 = 11.894 AU the origin is about to go stable again, (d) At 02 = 13 AU we are past the seeond 
bifureation; there is an inner eireulating zone surrounded by two libration lobes. 












































Extended Data Figure 5 
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Extended Data Figure 5: Characteristics of the LLER Islands: (a) Drift of EEER centers and 
extremities (again r] is rescaled by a factor of to turn it into an eccentricity-like 

variable), (h) The eccentricity of the outer planet increases significantly during capture in EEER. 
This is the dashed red line in Fig.l which is compared with the fiducial iV-wire simulation, (c) 
Modest growth of the eccentricity of the inner planet when captured in EEER. (d) The adiabaticity 
index T plotted versus Cmax for r = 10^ Tb. 









